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ABSTRACT 

The observed properties of galaxies vary with inchnation; for most apphcations we would rather have 
properties that are independent of inclination, intrinsic properties. One way to determine inclination 
corrections is to consider a large sample of galaxies, study how the observed properties of these 
galaxies depend on inclination and then remove this dependence to recover the intrinsic properties. 
We perform such an analysis for galaxies selected from the Sloan Digital Sky Survey which have been 
matched to galaxies from the Two-Micron All Sky Survey. We determine inclination corrections for 
these galaxies as a function of galaxy luminosity and Sersic index. In the g-band these corrections 
reach as as high as 1.2 mag and have a median value of 0.3 mag for all galaxies in our sample. We 
find that the corrections show little dependence on galaxy luminosity, except in the u band, but are 
strongly dependent on galaxy Sersic index. 

We find that the ratio of red-to-blue galaxies changes from 1:1 to 1:2 when going from observed 
to intrinsic colors for galaxies in the range —22.75 < Mk < —17.75. We also discuss how survey 
completeness and photometric redshifts should be determined when taking into account that observed 
and intrinsic properties differ. Finally, we examine whether previous determinations of stellar mass 
give an intrinsic quantity or one that depends on galaxy inclination. 

Subject headings: galaxies: clusters: general — galaxies: statistics — methods:statistical:surveys 



1. INTRODUCTION 

In our search to understand the formation and evo- 
lution of galaxies, some of our primary tools are mea- 
surements of the distribution of galaxy properties and 
relationships among these properties. From observations 
of these quantities at different redshifts we can deduce 
the nature of galactic evolution and by comparing them 
to the properties of dark matter halos we can constrain 
models of galaxy formation. 

The measurement of galactic distributions includes, 
but i s not limited to: the galaxy luminosity func- 
tion f Hubblc 193 6|), th e galaxy correlation function 
(jPeeble s & Hauscr, 11974 the distribution of galax- 
ies' s patial separatio i is), t he galaxy velocity func- 
tion (jGonzalez eta Ll l2000l), and the distributions 
of gala xy sizes ( Choloniewski 19 8j), s u rface bright- 
nesses (jFreemanl [T97a). colors (iBauml 119591: iFaljeij 
Il973f l. metalicit ies (lOsterbrockl 1 19701) and star for- 



mation r ates (jTinslev fc Danlvl ll980D . Relation- 
ships between ga laxy properties inclu de, the TuUy- 
Fisher re lation (iTullv fc Fish"ei1 Il977t ). the Faber- 
Jackson (Faber & Jackson' '1976") and fundamental 
plane (jDjorgovski & Davis 1987; Dressier ct al. 1987) 
relati ons, the luminosity-size relation (Kormendv 



1977f). the luminosity-metalicity relation p'aben il973i : 



Leaueux et al. 197^ and the density-morphology rela- 
tion (l Dresslei1ll980ii ). 

However, with the notable exception of the TuUy- 
Fisher relation these distributions and relations are tra- 



ditionally measured in terms of the o6ser?;e(i properties of 
galaxies. That is, the measurements used are K-corrected 
and corrected for foreground dust extinction, but no cor- 
rection is attempted to compensate for the viewing an- 
gle from which the galaxies are observed. In contrast, 
the TuUy-Fisher relation is not a relationship between a 
galaxy's observed luminosity and rotation velocity, but a 
relation between a galaxy luminosity and rotation veloc- 
ity corrected for inclination. The inclination correction 
attempts to recover the intrinsic properties of a galaxy 
and not properties that are measured because of the par- 
ticular angle from which the galaxy is viewed. Spiral 
galaxies are observed to have redder colors when their 
disks are more inclined, which is expected if the incli- 
nation increases the amount of dust that light traverses 
when emitted from the galaxy. 

Clearly, we would prefer to measure all galactic dis- 
tributions and relationships in terms of intrinsic galaxy 
properties instead of observed ones. The comparison of 
theory to observations is complicated and often done 
incorrectly because of confusion between observed and 
intrinsic galaxy properties. Early semi-analytic mod- 
els were unable to match both the galaxy luminosity 
function and the TuUy-Fisher relation in part because 
they failed to take into account that the first is ob- 
ser ved luminosity while the seco nd is intrinsic luminos- 
ity (jSomerville Sz Primacld[T999l ) . Also, when comparing 
galaxies at different redshifts we would like to be able to 
distinguish between evolution in their stellar populations 
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and changes in their dust properties. 

Furthermore, inchnation effects are of great help in un- 
derstanding the nature of dust in galaxies. Theoretical 
modeling of attenuation in galaxies is complicated be- 
cause it not only depends on the properties of dust, which 
seem to vary between galaxies, but also on how the dust 
is distributed and mixed with stars. By determining the 
intrinsic properties of galaxies we also learn how those 
properties change as a function of galaxy inclination and 
therefore some properties of the dust distribution. 

To determine intrinsic properties we need to know how 
a galaxy's properties change as a function of its inclina- 
tion. This is different then removing the effects of dust 
and dust is still present for a face-on galaxy. There are 
a number of approaches for addressing this issue each 
of which has its own merits and disadvantages. One 
approach is to solve for an inclination correction that 
minimizes the scatter in the TuUy-Fisher relation (e.g., 
lVerheiienll200lD . which assumes that the scatter in this 
relation should be as small as possible. Another method 
is to fit stellar population models to the SED of a galaxy 
and then assume that any discre pancies are caused by 
dust (e.g.. lKauffmann et al.ll2003h . A third is to observe 
background objects behind a foreground galaxy to get 
a direct measure of the extinction through the galaxy 
(e.g.,[Berhnd ct al. 1997; Holwerda et al. 2005), but this 
is difficult to do for more than a handful of cases. Fi- 
nally, one c an simulate the rad iative transfer through a 
galaxy (e.g.. lRocha et al.ll2007[ ). assuming one knows the 
distribution and scattering properties of the dust. 

The approach we explore here is somewhat simpler in 
that it assumes no knowledge of stellar population or 
dust properties. Instead, the main assumption is that 
a galaxy's properties should be independent of inclina- 
tion. Thus any statistical correlation between a galaxy 
property and inclination can be attributed to dust and 
the inclination correction is whatever makes the observed 
correlation go a way. This procedure has been applied a 
numb er of times (iGiovanelli et al.|[l99lll995t lT ullv et all 
119981 [Masters et all 120031 : IShao et al.ll200'7f K In this pa- 
per we greatly expand upon this method by applying 
it to 10,340 g alaxies taken from the Sloan Digital Sky 
Survey (SDSS, lYork et al.l [20001 with accompanying in- 
frared m agnitudes from the Tw o-Micron All Sky Survey 
f2MASS. ISkrutskie et "alll2006[ ). It is important to have 
galaxies with near infrared photometry because the ef- 
fects of attenuation a re minimized in these wavebands 
(jBell fc de Jong|[2001h . In a subsequent paper we will ex- 
tend the analysis performed here to the full SDSS galaxy 
catalog. 

We describe the method for determining inclination 
corrections in [2) In SJ3] we describe the sample we will 
use and discuss some of the properties of galaxies in this 
sample. In 2] we determine inclination corrections using 
our sample and compare our results to other determina- 
tions in fj5l In Sj6]we discuss how consideration of intrinsic 
properties can change our conclusions about the distribu- 
tion of galaxy properties focusing on the color- magnitude 
diagram. We also comment on the effect on survey com- 
pleteness, stellar masses and photometric redshifts. ij7] 
contains our conclusions and some discussion of future 
directions. 

2. THE METHOD 
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Fig. 1. — Histograms of visually classified elliptical (shaded) and 
disk (line) galaxies are shown in Sersic index, ris, g — r color and 
axis ratio, b/a. In each case the distribution of elliptical and disk 
galaxies are very different. Often Sersic index or color is used to 
divide elliptical and disk galaxies, but as is clear from the figure 
this still leaves many disks galaxies on the wrong side of the divide. 
Axis ratio can be used as a strong discriminator of galaxy type and, 
when combined with Sersic index, gives a sample that includes 70% 
of all disk galaxies. 

The method of determining extinction corrections sta- 
tistically is based on looking for correlations between a 
galaxy property and galaxy inclination. Under the as- 
sumption that the intrinsic properties of galaxies do not 
depend on inclination, one can infer the effect of at- 
tenuation by plotting how an observed galaxy property 
changes with inclination. This is shown in Figure [5] for 
the case of galaxy color, the property we will focus on in 
this paper. The extinction correction then is whatever is 
needed to remove the observed correlation between the 
chosen property and galaxy inclination. Note that this 
procedure only measures attenuation relative to face-on 
galaxies, it can not say anything about the total attenu- 
ation that occurs in a galaxy. ^ 

This method has been applied to late type galaxies 
where the property observed to vary with inclination is 
magnitude in an isophotal radius, galaxy color, or the 
galaxy luminosity function (IGiov anelli et al.lll994l. Il995t 
iTullv et all Hooa [Masters et al.l l2003: Sh ao et al.ll20?m . 
In order to describe these various studies within one 
framework we turn to a more mathematical description 
of the procedure. 

Galaxy surveys measure a number of galaxy proper- 
ties; fluxes, redshift, surface brightness, half light ra- 
dius, axis ratio, position angle, etc. Fluxes are usually 
converted into luminosities, a quantity that doesn't de- 
pend on the galaxy's distance (it has translational in- 
variance) using the galaxy's redshift and a K-correction. 
These constitute a set of observed galaxy data, G° = 
{M^,^°,r?^Q,b/a,pa, etc}. However, we would like to 
know the intrinsic properties of the galaxy, G*, those 
properties that are invariant to rotation and translation. 
We assume that the two are related by a transformation, 

^ It is also possible that an observed correlation between a 
galaxy property and inclination is the result of the data reduction. 
Nonetheless, in this case we would still like to correct this bias. 
The interpretation that the correlation is caused by dust therefore 
does require study by analyzing the data pipeline and comparing 
to expectations from dust modeling. For simplicity, in this paper 
we will assume that dust is the dominate source of any observed 
correlations. 
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T, such that 



(1) 



where T depends on the inchnation angle of the galaxy, 
9, and possibly on other galaxy properties. If intrinsic 
galaxy properties do not depend on inclination then we 
can determine the intrinsic galaxy properties by solving 
for the inverse transformation as the operation that 
satisfies, 



d{T-^G°) 



de 







G' 



(2) 



Note that this derivative is with respect to intrinsic 
galaxy properties. In practice this reduces to assuming 
a functional form for and then optimizing the pa- 
rameters in that function to come as close as possible to 
equation [2l 

One's choice about the functional form of T^^ there- 
fore ends up having a strong effect on the conclusions 
reached. Often it is assumed that the attenuation in a 
given waveband, A\, is of the form 

= -7Alog(cos0) (3) 

(jGiovanelli et al.lll994l . and others), where 7a represents 
the combination of dust properties, the distribution of 
dust and the distribution of stars that determines the 
attenuation at a given inclination. However, there is ev- 
idence that the dependence on inclination may be bet- 
ter fi t by log^(cos6') (jMasters et al.l l2003t iRocha et al] 
l2007f) . The parameter gamma was taken to be a constant 
for a given wavelength in early work (jGiovanelli et all 
|1994|) . but subseque ntly has been considered as a func- 
tion of luminosity (iGiovanelli et al.l 119951: iTuUv et ahl 
119981: [Masters et all 120031 : iShao et all 120071 ^. However, 
without a full understanding of the nature and distribu- 
tion of dust in galaxies we can hardly know what galaxy 
properties 7a should depend on. We must rely on the 
data and explore what functional forms and dependen- 
cies fit best. 

3. THE DATA 

Our g alaxy sample has be en taken from the NYU- 
VAGC (jBlanton et al.l l2003bD . We use galaxies from 
SDSS-DR2 that have near infrared magnitudes from 
2MASS (for a discussion of t he nature and co mplete- 
ness of these galaxies see ,McIntosh et al.l |2006[ ) . Each 
galaxy has been fit in the r-band with a Sersic profile 
using elliptical isophotes giving us a Sersic index, n^, a 
half light radius r^o and an axis-ratio b/a. We also have 
u, g, r, i and z total magnitudes from SDSS and J, H and 
Kg total magnitudes from 2MASS. In order to insure 
that our inclination measurements are accurate we re- 
strict our sample to galaxies with a seeing-deconvolved 
half light radius of three or more pixels. We also restrict 
ourselves to fits that are within the bounds of the allowed 
parameter space and not at one of the limits. Thus we 
only include galaxies with Sersic index and inclination in 
the range 0.5 < Us < 6.0 and 0.15 < b/a < 1.0. Finally 
we restrict our sample to —17.75 > Mk < —22.75 and 
'''50 £ 15kpc, the range where we have enough galaxies 
to make statistical statements about our sample. With 
these cuts we are left with 10, 340 galaxies. Furthermore, 
1,634 of these galaxies with Mr < —20.6 have been vi- 
sually classified by one of us (MRB). Because of the r 



TABLE 1 

Fit parameters to mean colors for face-on galaxies. 



mean color 



VQ 



VK 



Sigma 



SK 



1.67 
0.71 

0.30 
0.06 
-0.02 



u^^J.K) -0.01 



-0.27 
-0.23 
-0.20 
-0.18 
-0.17 
-0.01 
-0.02 



0.26 
0.10 
0.02 
0.01 
-0.05 
0.01 
0.00 



^(J-K) 



0.18 0.00 0.28 

0.25 0.07 0.08 

0.28 0.06 0.04 

0.28 0.06 0.02 

0.28 0.04 0.01 

0.12 0.04 0.01 

0.13 0.06 0.01 



Note. — This table lists the values of the parameters used to fit 
the mean distribution of faco-on galaxy colors according to equation 
[6]and the standard deviation about that mean according to equation 

m 

selection of SDSS and K selection of 2MASS this is not 
a fiux limited sample and thus we can not discuss the 
space densities of these galaxies. Therefore we will not 
discuss space volume at all in this paper, but will defer 
such discussion to a following paper that expands on our 
treatment here to the full SDSS galaxy catalog. 

3.1. Disk Galaxies 

As described in f2]the method is based on identifying 
observed properties that depend on inclination and then 
removing this dependency. However, inclination is not a 
directly observable galaxy property. Instead, what one 
measures is axis ratio. Previous studies have attempted 
to convert axis ratio to inclination assuming a thickness 
for a galaxy's disk. Before one can even do this though, 
one has to first identify which galaxies are disk galaxies. 
Usually this is done by cutting the sample by concen- 
tration, Sersic index or color. Fig [T] shows histograms of 
the 1634 visually classified galaxies with Mr < —20.6 in 
Sersic index, g — r color and axis ratio separately for disk 
(solid fine) and elliptical (shaded region) galaxies. We 
include SO galaxies as disk galaxies because they have a 
disk and therefore their measured axis ratio should be 
more related to the disks inclination then to the ellip- 
ticity of the spheroid component. It is clear that almost 
pure disk samples can be gotten by restricting one's sam- 
ple to Us < 3.0 or g — r < 0.7; however, in both these 
cases close to half of the disk galaxies are not included in 
the sample. Combining these two requirements leaves to 
a very mild improvement in sample completeness of only 
~ 5%. This is surprising since such cuts are often used 
in the literature to separate early and late type galaxy 
populations. It is worth noting that when such a cut is 
employed roughly half of the galaxies designated as early 
type in our sample arc in fact disk galaxies. 

Another property that shows a strong difference be- 
tween disk and elliptical galaxies is axis ratio. Elliptical 
galaxies almost never have small axis ratios and galax- 
ies with b/a < 0.55 are 90% disk galaxies in our sample. 
This is understandable as disks are intrinsically thin and 
thus can have very small axis ratios when seen in projec- 
tion, but ellipticals are close to spheroids whose intrinsic 
axis ratio qz, is rarely less than 0.5 and will be larger 
than this when seen in projection. Thus axis ratio can 
be used to determine if a red concentrated galaxy is truly 
an elliptical galaxy or not. Since red concentrated galax- 
ies are almost always assumed to be early-type in the 
literature one may be skeptical of our claim that many 
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Fig. 2. — Images of red {g — r > 0.75) and concentrated (ris > 4.0) and inclined (b/a < 0.5) galaxies from our sample. Clearly these red, 
concentrated galaxies are disk galaxies. Besides the cuts stated the only selection in these images is that the galaxy has a large angular 
size so that it makes a nice image. Axis ratio is an important diagnostic of galaxy type. 



of these galaxies are disk galaxies. Thus we show an 
example of eighteen such galaxies in Figure [2l These 
galaxies are selected to have ng > 4.0 and g — r > 0.75 
and b/a < 0.5. These are the eighteen largest galax- 
ies (in order to make nice postage stamps) that meet 
the above criteria. Clearly, highly-inclined, red, concen- 
trated galaxies are not elliptical galaxies. This is some- 
thing that should be taken into account when one is try- 
ing to identify elliptical galaxies. 

If we select all galaxies with Ug < 3.0 or b/a < 0.55 this 
gives a sample that is 94% disk galaxies and includes 70% 
of all disk galaxies. The disk galaxies missing from our 
sample are concentrated face-on disks. These galaxies 
will have the smallest inclination corrections and there- 
fore are of the least concern for our application. However, 
we need these galaxies to determine if the properties of 
concentrated disk galaxies are changing with inclination. 
We discuss how we deal with this issue in 21 

3.2. Inclination from axis ratio 

As mentioned above, the derivative in equation [2] is 
with respect to inclination, but inclination is not a di- 
rectly measurable quantity. Previous studies have dealt 
with this by assuming that all disks have some average 
thickness and then turning the measured axis ratio into 
an inclination by 

b/a = \/5.2 + (l-g2)cos(0), (4) 
where Qz is the ratio between vertical and radial scale 
heights. This equation should hold if the three- 
dimensional light distribution is well fit by concentric 
ellipsoids and the disk is optically thin. However, if one 
examines Figure [31 which shows the distribution of ob- 
served axis ratios for galaxies with < 3.0, binned by 
Sersic index, one sees that only the lowest Ug galaxies 
come close to having the distribution expected for ran- 
domly inclined disks with Qz = 0.15. All bins show a 



deficit of face-on galaxies, which is easily understood as 
all sources of asymmetry in the galaxy will push one away 
from perfectly circular isophotes. However, for galaxies 
with Us > 1.2 the existence of a bulge prevents these 
galaxies from having very low axis ratios. We do see 
that galaxies with Ug > 1.2 show essentially the same 
distribution of 6/a, implying that the measured b/a is 
mostly correlated with galaxy inclination. Because of 
these measurement issues we do not believe there is a 
reliable way to infer a galaxy's inclination from its axis 
ratio. It is possible that bulge-disk decomposition fits 
may yield better results for the disk inclination and thus 
reduce this source of uncertainty. For this paper we will 
simply work directly with axis ratio and reform equa- 
tion[2]to be that a galaxy's intrinsic properties shouldn't 
depend on observed axis ratio. 

4. DEPENDENCE OF OBSERVED PROPERTIES ON AXIS 

RATIO 

We now turn to studying the dependence of observed 
galaxy properties on axis ratio. We will start with galaxy 
color, which is the property that we have found to have 
the strongest dependence on axis ratio. One complica- 
tion when applying the method described in ij2] to the 
property of color is that since color is a difference be- 
tween two galaxy properties the correction found will not 
include any correction needed to the longer wavelength 
magnitude. As the K-hand is the longest wavelength we 
have access to we will consider X — K colors to minimize 
this effect. Thus what is recovered in our case will be 

= Ax^tot - Ak " Ax,e=o (5) 

where A\ tot is the total attenuation in a given waveband, 
Ak is any attenuation in the K-band and ^A,e=o is any 
attenuation in the face-on configuration. This is why it 
is important to have infrared magnitudes for our sample, 
as we expect Ak to be relatively small (|Bell fc de Jongj 



5 




0.02 - f - 

"i E 

0.00 L__L 

0.0 0.2 0.4 0.6 0.8 1.0 

b/a 

Fig. 3. — The distribution of axis ratios, for galaxies in our sam- 
ple in bins of Sersic index. The shaded region shows the theoreti- 
cal expectation for the axis ratio of a randomly inclined flattened 
spheroid with intrinsic axis ratio, = 0.15. For all bins there is 
a deficiency of nearly circular galaxies compared to theory as is 
expected since any asymmetry in the light profile will move the 
isophotes away from circular. Galaxies with rig > 1.2 have nearly 
the same distribution of axis ratios. Only the least concentrated 
galaxies have an axis ratio distribution similar to that of randomly 
oriented disks with finite thickness. Even a small bulge reduces the 
ellipticity of a galaxy's isophotes. Nonetheless, it is evident that 
l axia n atio is mpstly a. measure of jjalaxy inclination. 

120011) . In this section we wilronly snow ngures lor g — K 
colors, but our analysis is done for all six X — K colors we 
can produce. To see the analogous figures for other wave- 
bands please look at the supplemental online material or 
visit the website www.galaxystats.com/iirtriiisic7l 

In order to determine how a given galaxy's color de- 
pends on its inclination we need to know what that 
galaxy's color should be if we could see it face-on. Fig- 
ure m shows the distribution g — K colors of face-on 
disk galaxies {b/a > 0.85) as a function of if-band lu- 
minosity and Sersic index. The first two panels show 
all face-on galaxies with rig < 3.0 while the third panel 
shows the face-on galaxies with Ug > 3.0 that have been 
visually classified as disk galaxies by one of us (MRB). 
We only include the visually classified galaxies because 
we do not want to include elliptical galaxies, even though 
this means we are sampling a different luminosity range 
then the first two panels. We see that color depends 
on both Sersic index and luminosity which are not inde- 
pendent of one another. The relationship looks close to 
linear for absolute magnitude and linear in Sersic index 
to a value of = 4.0 at which point it appears to flatten 
out. Based on this we will assume that the dependence 
on Mk and rig are seperable and fit the mean galaxy 
color by a function of the form 

i^(A-/f) =Vo + Vk{Mk + 20) + VnUcS (6) 

where vi^x-K) is the average M\ — Mk color of a face-on 
galaxy with absolute magnitude Mk and Sersic index Ug 
and Tioff — rig for Ug < 4.0 and rice — 4.0 for > 4.0. 
The values of vq , vk and u„ that best fit the data are 
shown in Table [TJ One sees that the dependence of color 
on Us is relatively weak except for u — K colors, so this 
truncation for < 4.0 generally has very little effect. 

As can be seen in Figure 21 the dispersion about this 
mean also varies with Mk and n^, so we also fit the 
standard deviation with the formula 

<^{\-K) = So + sk{Mk + 20) + s„noff'. (7) 
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Fig. 4. — The g-K color of face-on galaxies {b/a > 0.85) is 
shown versus Mj^ in the left panel and ris in the right two panels. 
Diamonds show the mean value in a bin. For ris > 3.0 (the right 
panel) we only plot galaxies visually classfied as disks to avoid any 
confusion with ellipticals. The sample in the right panel is thus 
different then in the first two panels. In spite of this, it is clear 
that galaxy color is a function of luminosity and Sersic index. In 
both cases the mean relationship is close to linear, with evidence 
for a fiattening for Us > 4.0. 

Again we use a value of riea = 4.0 for Ug > 4.0 galaxies 
just to be consistent with equation [6l Values of the fit 
parameters for dispersion are also given in Table [TJ One 
can see that the dependence on rig is weak and whether 
or not one flattens the relationship at Ug — 4.0 will have 
very little effect on our results. In all cases, we find that 
these simple linear functions provides an adequate fit to 
the mean galaxy color and the standard deviation about 
this mean. 

With a model for the mean color of a face-on galaxy, 
we can now investigate how galaxy colors deviate from 
their mean face-on value as a function of axis ratio. Fig- 
ure [5] shows g — K color minus V(g-K) versus log6/a, 
for different bins of Sersic index. In all cases, a strong 
relationship is seen and one that seems to be roughly lin- 
ear in log6/a. The large intrinsic scatter in galaxy color 
makes it difficult to judge whether a more complicated 
dependence on log&/a is warranted. As a linear relation- 
ship seems adequate, we will express the attenuation in 
the form 

Ax = --ix\ogh/a. (8) 

Note that this is different than what has been done his- 
torically where the attenuation is considered to be linear 
in log (cos 0) (equation [3]). Here we are expressing the 
attenuation in terms of the observed axis ratio. It would 
be very interesting to compare this to the dependence of 
attenuation on observed axis ratio in radiative transfer 
simulations. 

Examining Figure [5l we see that 7^ (the slope of the 
line) varies with Sersic index. Creating plots like Figure 
Eland dividing the sample by other galaxy properties, we 
find that the correlation between color and axis ratio is 
strongest when binned by Ug or Mk- Therefore we un- 
cover that ^\ depends on at least these two parameters. 
We start with the simplest assumption, that 7^ can be 
determined from a linear combination of these parame- 
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Fig. 5. — The difference between a galaxy's g — K color and the 
mean value of the g — K color for a face-on galaxy with the same 
Mx and ris is shown as a function of axis ratio for six different bins 
of ris- We see in all cases, even for the highest n bin, there is a 
clear trend for more inclined galaxies to be redder. The relationship 
seems to be linear in logfe/a with a slope that becomes shallower 
with higher n. 



ters: 



7a = ao + aKiMx + 20) + a„ncff, 



(9) 



where the parameters ao,a„ and depend on wave- 
length. Again we will assume that dependence flattens 
off for Us > 4.0. Note that previous studies have only 
considered a linear dependence on luminosity, not Sersic- 
index. Now that we have a model for the attenuation we 
can investigate the best fit parameters for equationlHl To 
do this we use a Monte Carlo Markov Chain (MCMC) 
to minimize the statistic, 



X 



E 



{Mx -Ax- Mk) - v{x-K) 



(10) 



where Ax is given by equations [8] and [9l The resulting 
best fit parameters are shown in Tabled Figure [5] shows 
the effect of using these corrections to determine the in- 
trinsic colors of galaxies in our sample. Plotted is the 
intrinsic g-K color versus axis ratio for our disk galaxy 
sample We see that, unlike in Figure[5l there is not a dis- 
cernible correlation between galaxy color and axis ratio. 
Thus, it seems that our simple linear functions are suffi- 
cient to describe the behavior of the attenuation. When 
Figure [6] is made for the other 5 wavebands, a similar 
lack of correlation is seen. 

Table [2] also gives the mean value of 7a, the mean value 
of the attenuation and the maximum attenuation of any 
galaxy in our sample for each waveband. We see that 
in the g band the attenuation can be as large as 1.25 
magnitudes, but that the mean attenuation is 0.28 mag- 
nitudes. These values decrease substantially with longer 
wavebands reaching a mean attenuation of only a few 
hundredths of a magnitude for the near infrared bands. 
We expect any attenuation in the K-h&iyd to be less than 
this and thus negligible. 

As a semi-independent check of our results, we compare 
the distribution of Mg magnitudes and g — r colors for 
rig < 3.0 galaxies in our sample. We restrict ourselves to 
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Fig. 6. — The difference between a galaxy's intrinsic g — K 
color (corrected for inclination) and the mean value of the g — K 
color for a face-on galaxy with the same Mx and Us is shown 
as a function of axis ratio for six different bins of ns. We see that 
applying the inclination correction removes the dependence of color 
on inclination seen in FigO 

these low concentration galaxies so that we can compare 
face-on {b/a > 0.85) and edge-on {b/a < 0.30) disk galax- 
ies without having contamination from ellipticals. The 
distributions of these quantities for for face-on (shaded) , 
uncorrected edge-on (dashed line) and corrected edge-on 
(solid line) galaxies is shown in Figure [71 One sees that 
the distributions of face-on and edge-on galaxies is sig- 
nificantly different. The mean Mg magnitude and g — r 
color of face-on galaxies are —19.2 and 0.54, respectively; 
for edge-on galaxies the values are —18.6 and 0.71. This 
significant difference is dramatically reduced when con- 
sidering the corrected magnitudes of the edge-on galax- 
ies, which have mean values of —19.3 and 0.56, respec- 
tively. Examination of the histograms in Figure [7] shows 
that while corrected g — K colors have nearly the same 
mean as the face-on colors, there is some difference in the 
shape of the distribution. This seems to be an indication 
that an additional parameter or possibly a cross term 
is needed to fully describe the correction, though this is 
difficult to tell since there is a small bias introduced be- 
cause edge-on galaxies are included in the sample based 
on their observed and not their intrinsic magnitudes. Ei- 
ther way, these differences are relatively small and the 
overall success of our inclination corrections is evident. 

Having determined the correction for magnitudes we 
now turn to another galaxy property that shows inclina- 
tion dependence, galaxy size. Our measured half light 
radii r^o are only in the r-band, though it would be very 
informative for dust modeling to look at the radial depen- 
dence of attenuation in different wavebands. Examining 
the face-on sizes of galaxies in our sample we see that 
they dep end on K-band l uminosity (i.e., the Kormandy 
relation. lKormendvlll977t ). but show no dependence on 
sersic index. We find the mean face-on half light radii 
are well fit by 



logrso = 0.5 - 0.13(Mk + 20). 



(11) 



The distribution around this v alue is roughly lognor- 
mal as seen in p revious studies (jde Jong fc Lacevll2000t 
iShen et al.l [20031 so we will perform our correction in 
logrsQ. Comparing edge-on to face-on galaxies we find 
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TABLE 2 

Inclination correction parameters for different wave 

BANDS 



Band 


QO 


Or 


a„ 


max ■yx 


mean Ax 


max Ax 




1.79 


-0.43 


-0.28 


1.31 


0.39 


1.88 




1.38 


-0.19 


-0.22 


0.94 


0.28 


1.24 


a' 


1.02 


-0.09 


-0.12 


0.78 


0.23 


0.88 


Ai 


0.90 


-0.05 


-0.14 


0.60 


0.18 


0.70 


A, 


0.56 


0.01 


-0.01 


0.53 


0.15 


0.46 


Aj 


0.31 


-0.02 


-0.07 


0.17 


0.05 


0.24 


Ah 


0.15 


0.01 


-0.05 


0.05 


0.02 


0.11 



Note. — The table shows the parameters used to fit for ■yx to 
determine the attenuation in different bands according to eqns. 
|8]and[9] Also shown are the mean and maximum attenuation 
in each waveband and the maximum value of 7;^. 

that galaxies with < 2.0 have larger half light radii 
when they are inclined. For galaxies with Us > 2.0 we see 
no effect, probably because of the prominence of bulges 
in these galaxies. This increase of the sizes of inclined 
galaxies is what one would expect if there is more dust 
in the inner parts of galaxies causing greater attenuation 
in the center. However, as noted earlier, this may also 
be an artificate of the data pipeline used to determine 
galaxy size or a combination of both dust and the data 
pipeline. We parametrize the effect for galaxies by the 
formula 



log ''50 = logrgo +/3r(log&/a) 



(12) 



where r^Q and r^Q are the intrinsic and observed half light 
radii, respectively, and Pr is the strength of the effect in 
the r-band (analogous to 7^). Minimizing gives val- 
ues of /3r = 0.2. We find that we do not need luminosity 
or sersic index dependence to bring the edge-on galaxies 
into agreement with the face-on galaxy size distribution. 
This can be seen in the third panel of Figure [7] which 
shows the face-on and the corrected edge-on galaxy size 
distributions to be in good agreement. It will be inter- 
esting to continue this analysis for the full SDSS galaxy 
catalogue as we plan to in a following paper. Having de- 
termined corrections for observed magnitudes and sizes 
of galaxies in our sample, we now turn to comparing our 
results to previous determinations of the inclination cor- 
rection in the literature. 

5. COMPARISON TO PREVIOUS WORK 

In comparing our results to previous work it is impor- 
tant to take into account three differences in our treat- 
ment. One, our value of 7^ is multiplied by log6/a in- 
stead of log (cos 6*). Two, we explore 7a dependence on 
both Mk and rig, while earlier studies only considered 
dependence on luminosity. Finally, we include all types 
of galaxies in our sample instead of focusing on just late 
type galaxies. 

Our values of 7^ for an = 1.0 galaxy are shown as a 
function of galaxy luminosity in Figure [5] One sees that 
aside from in the u-band, the dependence on luminos- 
ity is fairly we ak. For comparison, we show the values 
dete rmined bv iTullv et al.l (|l998f ): [Masters et all (|2003l ) 



and lShao et al.l (|2007D ." Note that IShao et al.l found no 
evidence for luminosity dependence, so thier fit assumes 
this and that 7a has a power law depend ence on wave- 
length. Also, to compare to iTuUv et al.l . we have con- 
verted the measured dependence to Mk by subtracting 




Fig. 7. — The left most panel shows the distribution of observed 
g-band magnitudes for face-on (shaded) and edge-on (dashed line) 
galaxies. Also shown is the distribution of intrinsic (corrected) 
magnitudes for the edge-on (solid line) galaxies. One sees that the 
distributions of observed magnitudes is very different, but that once 
an inclination correction has been applied the distributions become 
almost identical. The middle panel shows distributions of g — r 
color, with the same line-styles as the left most panel. Again the 
inclination corrections bring the edge-on galaxies into agreement 
with the face-on ones. However, the shape of the corrected color 
distribution differs from the face-on one, suggesting that their may 
be dependence on another parameter needed for the correction. 
The right most panel shows the distribution of observed galaxy half 
light radii, rso, again with the same line styles indicating observed 
face-on, observed edge-on and intrinsic edge-on galaxies. In all 
three panels only galaxies with < 3.0 are used to insure a fair 
comparison. 

the mean r — K and i — K colors of our galaxy sample. 
We see that the determinations of 7a from all studies 
have similar values, that is they all agree with one an- 
other at some luminosity. However, the dependence on 
luminosity shows a wide range of behaviors. There are 
many possible reasons for these differences. One is defi- 
nitely the assumpt ions that have been made, for example 
the 77 value from [Masters et al.l assumes no luminosity 
dependence, while the 7j values is fit to a broken line. 
Clearly, we expect that we will get different results be- 
cause we have fit to both luminosity and Sersic index. 
Other issues that may be important are the sample se- 
lection used in the study, the conversion from b/a to 
inclination done in the other studies, a nd which p rop- 
erty is being correlated with inclination. 'TuUv et al.l and 
we both used galaxy color, while Masters et a L looked at 
magnitude in an isophotal radius and lShao et al.l focused 
on the galaxy luminosity function. In particular, as we 
have shown in SJU not only is the total luminosity of a 
galaxy dependent on axis ratio, but also the half light 
radius is too. Since isophotal radius depends on both it 
may not be surprisi ng that we arrive at different results 
than iMasters et al.l These intricacies will require further 
study, which would be most enlightening when combined 
with radiative transfer simulations. 

The dependence of 7a on Sersic index is shown in Fig- 
ure [9] In general, we find that the dependence on Sersic 
index is stronger than on K-hand magnitude. We have 
fit each band independently so the unphysical situation 
that arises in Figure [5] can occur where 7^ becomes larger 
than 7i at high Sersic index. This is of course undesir- 
able, but our approach is to investigate what we can learn 
from the data instead of coming in with theoretical preju- 
dices. One can easily find an acceptable solution to this 
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Fig. 8. — The value of ■yx for spiral galaxies (n^ = 1) as 
a function o f Mk is shown in comparison to other results from 
ITuUv et al] iITQQSI. da shed lines). Masters et al. (2003, dash-dot 
lines') and IShao et al.l 1 120011 . dotted lines). Note that none of the 
wavebands are the same (except for the J-band which are all from 
2MASS) and that the methods used to determine -yx all differ. 

by requiring that 7a is always smaller in longer wave- 
lengths. That would mean that in Figure [HI would 
take on the value of where now jz > 7i- One way to 
improve in this area would be to fit multiple wavelengths 
at the same time, under some assumptions about how 
the attenuation at different wavelengths is related. 

Since we fit each waveband independently and since 7a 
depends on both Mk and Us , the wavelength dependence 
on 7a, and thus the attenuation, varies with Mk and 
Hg. Figure [TOl shows the wavelength dependence of 7a 
for fixed values of Mk and n^. We see that generally a 
power law is a good fit to 7a in the SDSS wavebands. 
The index of the power law seems to increase with Mk 
for fixed ns, while at fixed ris the slope of the index of 
the power law remains almost unchanged and only the 
amplitude of 7a increases with decreasing n^. The power 
law index ranges from ~ 0.5 to ~ 1.5. These relationships 
should be very useful in understanding the distribution of 
dust in these galaxies. We see in all cases 7a in the near 
infrared bands (J and H) falls below the power law fit. It 
is clear from this that the amount of attenuation in the 
K-band, not included in our calculations, is negligible. 

We conclude from this that the corrections we find are 
in general agreement with previous results and that dif- 
ferences are probably primarily due to our including the 
dependence on Sersic index. We now turn to a prelim- 
inary exploration of how galaxy properties differ when 
considering the intrinsic instead of observed values. 

6. INTRINSIC GALAXY PROPERTIES COMPARED TO 
OBSERVED ONES 

In this section we turn to comparing intrinsic galaxy 
properties to observed ones. A complete discussion of 
this subject is beyond the scope of this work. And as 
we have mentioned earlier, we will not be discussing the 
space density of objects because we do not have a flux 
limited sample. Instead we will focus on how galaxy 
photometric quantities change in our sample when going 
from observed to intrinsic properties. The overall effect 
is fairly modest, as can be seen from Table [21 the aver- 
age attenuation in the g-band is Ag — 0.28. This will 



Fig. 9. — The value of 7^ as a function of Sersic index, -yx 
shows a strong dependence on Sersic index in all bands except for 
the z-band. 7^ shows almost no dependence on ris, which becomes 
unphysical when -jz > -fi- This seems to imply that our assumption 
that 7;^ is linear in Us should be modified. Unfortunately, unlike 
in Fig. Owe can not compare our results as other studies have not 
calculated -yx as a function of Sersic index. 

have small effects on the galaxy luminosity function and 
correlation function, essentially a small shift in their am- 
plitudes. However, the inclination correction is not uni- 
form, and while many galaxies have no correction others 
are brightened by up to 1.25 magnitudes in the f;-band. 
Thus the most notable effect will be in comparing the 
properties of spiral and elliptical galaxies. Wc will focus 
on the differences between observed and intrinsic galaxy 
color as seen in the color magnitude diagram, the de- 
termination of stellar masses, and the determination of 
photometric redshifts, all measurements that are based 
on galaxy color. First though, we start with discussing 
survey completeness. 

6.1. survey completeness 

When information is presented from a galaxy survey it 
is usually taken from a flux limited catalog. The SDSS 
spectroscopic sample is complete to an apparent magni- 
tude of rur = 17.7. This can be converted to an equiv- 
alent completeness in observed absolute magnitude at 
some distance. But this does not translate into an equiv- 
alent completeness in intrinsic magnitudes. Thus, when 
one looks at a flux limited catalog, there are galaxies 
that are intrinsically identical to galaxies in the catalog, 
but are excluded because of the inclination they hap- 
pen to be viewed at. When discussing the distribution 
of galaxy properties or connecting galaxies to dark mat- 
ter halos this is not desirable. In intrinsic magnitudes, 
the survey is only complete to the maximum attenuation 
brighter than the observed magnitude. So in the r-band, 
where the maximum attenuation is ~ 0.9 magnitudes, 
this corresponds to rrir — 16.8 a much shallower survey. 

If one wants an intrinsically complete survey it is not 
necessary to obtain spectra for every galaxy fainter than 
the survey limit down to this maximum attenuation. In- 
stead one only needs to target galaxies where 

< 17.7 -I- max(7r) * log &/a (13) 

a much smaller sample of galaxies. If one determines 
Sersic indeces before targeting the galaxies then max(7r) 
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Fig. 10. — The value of 7;^ as a function of wavelength for various 
values of Af/f at fixed rts (top panel) and for various values of Us 
at fixed Mfc (bottom panel). The lines are power law fits to the 5 
values for the SDSS bands. One sees that though the first 5 bands 
seem to be well fit by a power law the measured values of ■yj and 
-yn are always well below the extrapolation of this eurve. 

would equal the maximum value of 7^ for that value of 
rifi. Without the Sersic indeces (or some equivalent mea- 
sure) one would have to use the maximum value of 7^ 
which for our sample is 1.1. To give an example, we use 
the DR4 sample of SDSS, which does not have determi- 
nations of Sersic indeces, but has measured axis ratios, 
though different then the ones we have been using in 
this paper (see These axis ratios are measured with 
exponential fits (labeled ab_exp in the catalog), but are 
reasonably close to the ones we have been using which 
are based on elliptical Sersic fits. In this sample there are 
92,752 galaxies with < 16.8. Going 0.9 mag deeper 
to rrir < 17.7 increases the sample by 242,103 galax- 
ies. However, if the axis ratio of the galaxy is considered 
then there are only 89,855 galaxies that may be brighter 
than 16.8 mags after correcting for inclination. Thus 
taking the axis ratio into account reduces the required 
follow up spectra by 63%. Future surveys may want to 
consider axis-ratio along with apparent magnitude when 
determine their galaxy sample. In this way, samples lim- 
ited to some intrinsic magnitude can be compiled which 
will greatly facilitate comparisons to theory and other 
observations. 

6.2. color-magnitude relationship 

One area where we see a significant difference be- 
tween observed and intrinsic properties is in the color- 
magnitude diagram. Figure fTTl shows a color- magnitude 
diagram for observed and intrinsic magnitudes for our 
sample. There is no l/Vmax weighting as explained ear- 
lier, which would increase the difference between the two 
plots. Clearly, the distribution of points has changed in 
the two panels. In the observed color-magnitude dia- 
gram 46% of galaxies have g — r > 0.7; for the intrinsic 
color-magnitude diagram this decreases to 32%. There- 
fore, in observed colors the split between red and blue is 
roughly 1:1, in intrinsic color this changes to 1:2. Asser- 
tions about how the galaxy population is evolving based 
on the number o f red or blue galaxies (|Faber et al.,.2005, : 
iBell et al]|2007l) should probably look at changes in in- 




FlG. 11. — The observed eolor-magnitude diagram (left panel) 
and the intrinsic color-magnitude diagram (right panel). In the 
observed color magnitude diagram there are many galaxies redder 
than the red sequence. These galaxies are mostly removed in the 
intrinsic eolor-magnitude diagram. Face-on (b/a > 0.85) and edge- 
on {b/a < 0.35) galaxies are color coded by purple and orange, 
respectively. In the left panel they occupy a completely different 
region of the diagram, but in the right panel they have similar 
distributions. 

trinsic instead of observed properties if they want to be 
sure that reddening isn't a large part of the effect they 
are measuring. 

6.3. stellar masses 

One intrinsic galaxy property that we have not dis- 
cussed yet is stellar mass. Stellar mass is not directly 
observable, but there are a number of methods that 
try to estimate the stellar mass based on spectroscopy 
or photometry, (|Bell_| £iD |2003t iBlanton et all l2003al : 
iKauffman n et al.ll2003t ). Clearly, stellar mass IS an in- 
trinsic property - the mass of a galaxy shouldn't change 
under rotation. To check that this is the case then, one 
would like to see that equation [2] holds for stellar mass. 
Figure[l2]shows histograms of the stellar mass for face-on 
[b/a > 0.85, shaded) and edge- on (b/a < 0.3, line ) galax- 
ies using the mass esti mates of IKauffmann et al.l ([2003), 
IBlanton et all (|2003aD and IBell et al.l (120 031. We see 
that in the estimates of IKauffmann et alT edge-on galax- 
ies are slightly less massive then face-on galaxies, that for 
IBlanton et al.l edge-on galaxies are slightly rnore massive 
then face-on galaxies and that for iBell et alj there is re- 
markably good agreement between the two . We should 
note that the method of IKauffmann et al] is based on 
spectra which, for the low redshift sample we are con- 
sidering will only probe the inner parts of a galaxy and 
therefore is more likely to be in error. Also, we are only 
testing that the determined stellar mass does not depend 
on inclination not that it is the correct value of the stellar 
mass. The point of this exercise is just to demonstrate 
that one should check that intrinsic quantities actually 
don't show dependence on inclination. 

6.4. photometric redshifts 

Finally, we turn to the issue of det ermining photo- 
metric redshifts (e.g.. iHogg et al.l[l998l ) for galaxies, or 
clusters found by id entifying red-sequence galaxies (e.g 
iKoester et al.l l2007f) , which is a subset of photometric 
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Fig. 12. — The distribution of stellar masses in our sample for facc-on (shaded) and edge-on (line) galaxies with Us < 3.0. Th e left panel 
shows stellar masses determined bylKauffmann et al. (2003), the middle panel by Blan ton et al . (20033-) and the right panel bv lBell et al.l 
112003 ). We see that in the left panel the edge-on galaxes are not corrected enough, in the middle panel they are corrected too much and 
for the right panel the correction is just right. Note that this only tests that there is not inclination dependence in the determination, not 
that the stellar masses are correct. 



redshifts. In recent years algorithms to determine a 
galaxy's redshift only using photometric redshifts have 
become very successful. These methods rely on the fact 
that there is a limited family of intrinsic SEDs that galax- 
ies have so that with enough wavelength coverage it is 
possible to determine both the galaxy's SED and red- 
shift with only photometric data. However, as we have 
demonstrated in this paper, the observed SED of a galaxy 
changes with inclination. It may still be possible to dis- 
entangle the intrinsic SED, reddening and redshift with 
observations in enough wavebands, but in cases where an 
axis ratio can be measured, photometeric redshifts ignore 
useful information. Thus we strongly recommend that in 
future determinations of photometric redshifts axis ratio 
be used to help reduce uncertainties. 

7. CONCLUSION 

This paper has focused on the intrinsic properties of 
galaxies, whi ch are routinely dete rmined for TuUy-Fisher 
studies (^e.g.. lPizagno et aLll2007i ). but not for other sta- 
tistical studies of galaxies. Intrinsic properties are invari- 
ant under changes in viewing angle, unlike observed prop- 
erties which would change if we could view a galaxy from 
a different vantage point. Our main goal in this paper 
has been to clarify the difference between observed and 
intrinsic properties and to suggest how observed proper- 
ties can be converted into intrinsic ones. 

The method we use in this paper is to identify an ob- 
served galaxy property that shows a correlation with axis 
ratio and then apply the necessary correction to remove 
this correlation. We find that both color and size show 



correlations with axis ratio. We are able to remove these 
correlations using simple linear formula that depend on 
i^T-band magnitude and Sersic index. We therefore can 
construct a galaxy catalog with intrinsic value for galaxy 
size and magnitude. 

There are many distributions and relationships that 
should be reconsidered in terms of the intrinsic prop- 
erties instead of the observed ones. We can not cover 
all of these in a paper of this length but we highlight 
a few points to suggest how things may change. We 
focus on the color magnitude diagram as an example. 
The observed color magnitude diagram shows a number 
of galaxies redder than the red sequence and a lack of 
bright blue galaxies. When we plot the intrinsic color 
magnitude diagram these ultra red galaxies are mostly 
removed and there are many more bright blue galaxies. 
We find that the ratio of blue to red galaxies changes 
from 1:1 to 2:1 for galaxies with absolute luminosities 
—23.75 > Mji > —17.75, a significant change. 

In a following paper, we apply the insights we have 
gained here to producing inclination corrections for the 
full SDSS catalog. Having intrinsic properties for this 
flux limited catalog will then allow us to address how 
volume densities are effected when going from observed 
to intrinsic quantities. 
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